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I.  COMMENTS  ON  POWER  VS.  CURRENT  CURVES 


In  the  single-mode  regime,  the  intracavity  intensity  of  the  single-mode  determines  the 
strength  of  the  parametric  interaction  with  the  sidebands.  Therefore,  we  would  like  to 
calculate  the  intracavity  intensity  from  the  measured  output  power.  We  remain  true  to  the 
distributed  loss  approximation,  for  which  the  output  power  is  given  by 
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where  am  =  ln[l/(i?ii?2)]/(2L),  the  length,  width,  and  height  of  the  cavity  are  L,  w,  and 
h,  and  the  time-average  intensity  of  the  single-mode  is  ( E 2)  =  2|£0|2.  We  are  assuming  a 
uniform  field  intensity  in  the  transverse  dimensions,  and  therefore  not  worrying  about  the 
transverse  overlap  factor.  We  can  rearrange  this  equation  for  the  intracavity  intensity 
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With  this  equation,  we  can  convert  the  measured  total  output  power  of  each  laser  into  the 
intracavity  intensity,  using  our  measured  values  of  the  refractive  index  neg  and  the  dephasing 
time  T2,  our  best  estimates  for  d  and  Ti,  and  in  the  case  of  the  HR/AR  laser  we  have  used 
Ri  =  1,  i?2  =  0.01.  The  result  is  plotted  in  Fig.  SI  as  a  function  of  J/Jth,  which  is  the  same 
as  the  inset  in  Fig.  1  of  the  main  text. 

The  theoretical  formula  for  the  intracavity  intensity  is 

p  —  1 


w  = 


(3) 
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where  p  =  weq/wth  is  the  pump  parameter.  We  emphasize  that  p  is  not  the  same  as  -J/Jth- 
The  slope  of  \£q\ 2  vs.  p  is  always  between  2/3  and  1,  depending  on  the  diffusion  parameter 
7 £).  The  reference  line  in  Fig.  SI  is  drawn  with  a  slope  of  one  to  indicate  that  each  of  the 
| So | 2  vs.  J/Jth  curves  has  a  slope  greater  than  one.  Therefore,  we  conclude  that  J/Jth  must 
underestimate  p.  One  factor  that  contributes  to  this  underestimation  is  the  transparency 
current  Jtrans:  a  fixed  amount  of  current  that  must  be  delivered  to  the  active  region  simply 
to  raise  the  inversion  from  a  negative  number  to  zero.  To  understand  this  simply,  suppose 
that  the  equilibrium  inversion  scales  like  weq  oc  J  —  Jtraris ,  and  that  Jtrans  remains  a  constant 
number  at  threshold  and  above.  Then  the  pump  parameter  p  =  weq/wth  is  expressed  in 
terms  of  J  as 

J  J  trans  /  ,  > 
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FIG.  SI.  The  measured  output  power  of  each  QCL  (from  both  facets)  is  converted  to  the  intracavity 
intensity  |£o|2,  and  plotted  against  J/Jth-  Each  curve  is  color-coded  to  indicate  the  range  over  which 
the  laser  operates  in  a  single-mode,  harmonic  state,  or  dense  state.  Note  that  the  quantity  |£o|2 
is  only  meaningful  in  the  single-mode  regime,  because  we  are  interested  in  the  intensity  of  the 
single-mode  before  the  harmonic  regime  sets  in.  A  line  with  a  slope  of  one  is  drawn  as  a  reference. 

For  example,  suppose  that  for  a  laser  with  Jt h  =  500  mA  the  harmonic  state  kicks  in  at 
550  mA,  or  J/Jth  =  1.1.  If  the  transparency  current  was  Jtrans  =  250  mA,  (in  other  words, 
half  of  the  threshold  current,  which  is  reasonable  for  QCLs),  then  the  pump  parameter 
at  the  harmonic  state  onset  would  be  p  =  (550  —  250)/(500  —  250)  =  1.2.  Thus,  J/Jth 
underestimates  p. 

A  more  rigorous  study  is  required  to  determine  Jtrans  for  each  laser,  which  can  be  done  by 
measuring  many  lasers  of  the  same  active  region  but  different  lengths.  Once  Jtrans  is  known, 
the  slope  of  |£o|2  vs.  p  should  fall  between  2/3  and  1  and  in  principle  a  value  for  yy;  can  be 
extracted,  allowing  one  to  quantify  the  amount  of  diffusion  present. 

II.  HYSTERESIS  OF  IV  CURVE 

The  IV  curve  of  device  DS-3.8  shown  in  Fig.  S2  demonstrates  a  hysteresis.  When  starting 
below  threshold  and  increasing  the  current  (red),  the  voltage  of  the  laser  decreases  (negative 
differential  resistance)  when  the  noisy  harmonic  state  transitions  to  the  dense  state  at  523 
mA.  (In  the  spectra  shown  in  Fig.  3(c)  of  the  manuscript,  this  transition  occurs  at  502  mA. 
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Current  (A) 

FIG.  S2.  The  IV  curve  of  DS-3.8  exhibits  a  hysteresis  as  the  current  is  increased  (red)  and  decreased 
(blue).  The  hysteresis  is  correlated  with  the  transition  from  the  noisy  harmonic  state  to  the  dense 
state. 


The  exact  current  at  which  the  transition  occurs  is  not  identical  from  one  experiment  to 
another,  but  occurs  predictably  within  a  range  of  about  20  mA.)  Once  the  laser  reaches  the 
dense  state  and  the  current  is  subsequently  decreased  (blue),  the  laser  persist  in  the  dense 
state  at  currents  below  523  mA.  This  is  responsible  for  the  hysteresis  loop.  For  the  same 
current,  the  voltage  is  4.8  mV  smaller  when  in  the  dense  state  than  in  the  noisy  harmonic 
state.  At  494  mA,  the  dense  state  transitions  to  the  noisy  harmonic  state,  and  the  two 
voltage  curves  overlap  again. 

The  lower  voltage  of  the  dense  state  indicates  a  larger  radiative  photocurrent,  which 
implies  that  the  output  power  of  the  laser  is  slightly  larger  in  the  dense  state  than  in  the 
harmonic  state,  for  the  same  current  pumping.  (While  we  could  have  measured  the  output 
power  to  demonstrate  this,  the  IV  measurement  is  more  sensitive.)  Thus,  the  dense  state  is 
more  efficient  at  extracting  photons  from  upper  state  electrons  and  is,  in  one  sense,  more 
stable.  Therefore,  once  the  laser  enters  the  dense  state,  it  likes  to  remain  there  even  as  the 
current  is  decreased. 
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(a)  LL-9.8:  increasing  current 


(d)  LL-9.8:  decreasing  current 
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(b)  TL-4.6:  increasing  current 
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(c)  DS-3.8:  increasing  current 
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(e)  TL-4.6:  decreasing  current 
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(f)  DS-3.8:  decreasing  current 
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FIG.  S3.  The  spectra  of  the  three  uncoated  devices  as  the  current  is  ramped  up  and  as  it  is  ramped 
down.  Of  the  six  plots,  only  (e)  and  (f)  were  not  included  in  the  main  text. 
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(a)  TL-4.6:HR/AR:  increasing  current 

1220  mA  , 

:1165  mA 

1160  mA 

-  i  1 

’"""  -jlkr 

JL  A 

|*t.  1.68 
„  1.60 

1.60: 

1050  mA 

1 

1  -  . 

1.44 

1000  mA 

- . , . H 

i  . 

1.38- 

950  mA 

::::T±= 

1.31 

900  mA 

,  ,  L  Jw  JL 

. 1 . r  . 

1.24 

885  mA 

JL  Jl-,46  fsr 

1.22_ 

880  mA 

. r: .  * . . l  . 

1.21: 

875  mA 

1.20_ 

735  mA 

j/Jtf 

i  =  1.01 

> 

io 

m 

■o 

o 


1220  mA 
1100  mA 
1070  mA 


1040  mA 
960  mA 
935  mA 
890  mA 
885  mA 


880  mA 
875mA 
735  mA 


(b)  TL-4.6:HR/AR:  decreasing  current 
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FIG.  S4.  The  spectra  of  TL-4.6:HR/AR  as  the  current  is  (a)  ramped  up  and  (b)  ramped  down. 
Only  plot  (b)  was  not  presented  in  the  main  text. 


III.  ADDITIONAL  DATA  ON  SPECTRAL  HYSTERESIS 


In  Fig.  S3,  the  spectra  of  the  three  uncoated  devices  are  plotted  as  the  current  is  increased 
starting  from  below  threshold,  and  also  as  the  current  is  decreased.  For  devices  LL-9.8  and 
TL-4.6,  once  the  dense  state  is  reached  it  persists  as  the  current  is  decreased  until  the  single¬ 
mode  state  reappears.  Device  DS-3.8  also  exhibits  noisy  harmonic  states  as  the  current  is 
decreased.  For  all  three  devices,  it  is  true  that  the  clean  harmonic  state  never  reappears 
once  the  laser  has  reached  the  dense  state. 

In  Fig.  S4,  the  spectra  of  TL-4.6:HR/AR  are  plotted  as  the  current  is  increased  and  then 
decreased.  As  the  current  is  decreased,  the  dense  state  persists  until  a  noisy  harmonic  state 
appears  at  1.29Jth.  Interestingly,  in  this  device-  unlike  the  uncoated  ones-the  clean  harmonic 
state  with  one  pair  of  sidebands  reappears  at  885  mA,  which  is  quite  close  to  the  instability 
threshold  of  880  mA  found  when  the  current  is  ramped  upwards.  The  reappearance  of  the 
harmonic  state  may  also  suggest  that  the  sidebands  are  in  the  parametric  enhancement 
regime,  in  addition  to  the  two  pieces  of  evidence  discussed  in  the  main  text;  namely,  the 
large  sideband  spacing  (46  FSR)  and  the  large  range  of  intracavity  power  over  which  the 
harmonic  state  persists. 
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IV.  THEORY:  SINGLE-MODE  SOLUTION 


This  section  gives  a  more  detailed  derivation  of  the  single-mode  solution,  including  the 
intracavity  power  as  a  function  of  pumping,  and  the  population  inversion  as  a  function  of 
position  and  pumping. 

For  a  two  level  system  with  upper  state  |a)  and  lower  state  b) ,  the  material  equations  in 
the  non-rotating  frame  and  the  field  equation  are 


dpab 
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dt 

d2E  1  d2E 
dz2  c2  dt 2 
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We  emphasize  that  these  equations  are  in  the  non-rotating  frame,  whereas  the  equations  we 
have  used  in  the  main  text  [1]  were  already  in  the  rotating  frame  and  the  RWA  had  already 
been  applied.  However,  since  we  are  here  dealing  with  two  counter-propagating  waves,  we 
chose  to  more  closely  follow  the  approach  in  [2],  We  make  the  following  ansatzes: 


E{z,  t)  =  4=  [£r{z,  +  sL(z,  t)e-dMz)  +  c.c]  (8) 

v  2 

Pab(z ,  t)  =  7}*r(z,  t)e~dut~kd  +  v*L(z,  t)e-^t+k^  (9) 

w(z,t )  =  w-oc{zi  t)  +  w2{z,  t)el2kz  +  W2(z,t)e~l2kz .  (10) 


(We  use  the  subscript  “DC”  rather  than  “0”  for  the  spatial  average  of  the  population  inver¬ 
sion,  wdc,  because  the  subscript  0  is  used  throughout  the  text  to  refer  to  the  primary  mode. 
No  such  ambiguity  occurs  for  the  subscript  “2.”)  Plugging  the  ansatzes  into  the  differen¬ 
tial  equations,  and  making  the  RWA  as  well  as  the  slowly-varying  envelope  approximation 
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where  k  =  2d/h,  a  =  NujTp-P  \J p/e /h  is  the  Beer  absorption  coefficient  of  the  material,  and 
A  =  ujba  —  ijj  is  the  detuning  of  the  held  from  the  atomic  resonance  frequency.  The  loss  term 
4  has  been  added  to  the  held  equation  heuristically,  and  in  this  context  it  represents  only 
the  waveguide  loss. 

We  solve  for  the  single-mode  solution  by  setting  the  time-derivatives  to  zero  and  the 
slowly- varying  envelope  functions  to  be  constants.  In  doing  so,  we  are  now  making  the 
distributed  loss  approximation  because  we  are  not  allowing  the  helds  to  grow  in  space. 
Thus,  4  must  now  be  taken  to  be  the  total  loss,  waveguide  plus  mirror  loss.  We  take  A  =  0 
for  simplicity,  because  the  single-mode  will  lase  very  close  to  the  peak  of  the  gain  spectrum. 
We  denote  the  steady-state  held  amplitudes  by  Sr  =  £l  =  So  and  hnd  the  LI  curve 

p  —  1 
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where  7d  =  (1  +  Ak2DTi)~l  is  the  diffusion  parameter.  The  steady-state  population  w0(z) 
is  given  by 
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V.  THEORY:  POPULATION  PULSATIONS 


This  section  gives  a  more  detailed  derivation  of  the  population  pulsations,  and  demon¬ 
strates  how  to  include  nonzero  detuning  A  and  GVD  into  the  formalism. 

We  begin  by  imagining  a  small  volume  of  dipoles  subject  to  a  spatially  uniform  Li-held 
to  develop  an  understanding  of  the  non-linear  effects  caused  by  the  Bloch  dynamics.  The 


electric  held  is  given  by 

E(t)  =  £(t)ei0Jt  +  c.c. 

(19) 

The  Bloch  equations  in 

the  rotating  wave  approximation  are 
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•  /  c*  c  W  We(l 

W  =  IK{C  (7-6(7  ) - — - 1 

T\ 

(21) 

where  a  is  the  off-diagonal  element  of  the  density  matrix  in  the  rotating  frame,  w  is  the 
population  inversion  (positive  when  inverted),  A  =  uiM  —  uj  is  the  detuning  between  the 
applied  field  and  the  resonant  frequency  of  the  two-level  system,  T\  is  the  (longitudinal) 
population  relaxation  time,  T2  is  the  (transverse)  dephasing  time,  k  =  2 d/h  is  the  coupling 
constant  where  d  is  the  dipole  matrix  element  (assumed  to  be  real)  and  h  is  Planck’s  constant, 
and  weq  is  the  equilibrium  population  inversion  in  the  absence  of  any  electric  held  which  is 
determined  by  the  pumping.  (Note  that  these  equations  are  identical  to  Eqs.  3. 19(a)- (c) 
in  [1],  except  that  we  have  allowed  £  to  be  complex  and  left  the  off-diagonal  component  of 
the  density  matrix  in  complex  notation  rather  than  writing  a  =  (u  +  iv)/2.)  With  these 
conventions,  the  macroscopic  polarization  P  (dipole  moment  per  volume)  in  a  region  with 
a  volume  density  of  N  dipoles  is  given  by 

P{t)  =  Ndaeiut  +  c.c.  (22) 


First,  we  consider  the  effect  of  a  monochromatic  held  at  frequency  u>,  obtained  from  Eqs. 
20-21  by  setting  £{t)  =  £0  and  all  time  derivatives  to  zero.  The  result  is  a  steady-state 
polarization  <j0  and  population  inversion  uy  given  by 

inT2 


(J0  = 


2(1  -  iAT2) 


w0£0 


w, 


w0  = 


eq 


1  + 


K2Tir2|£o 


(23) 

(24) 


l+(AT2)2 

Note  that  the  population  inversion  uy  is  saturated  as  the  held  strength  £0  increases:  this  is 
responsible  for  saturable  loss  (when  weq  <  0)  and  saturable  gain  (when  weq  >  0). 


A.  Two-frequency  operation 

Next,  we  consider  the  E-field 

8  =  So  +  £+ei5ut  (25) 
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which  consists  of  the  strong  field  So  at  frequency  u>  superposed  with  the  much  weaker  field 
£+  detuned  from  uj  by  5uj.  A  polarization  will  of  course  be  induced  at  uj  +  5uj.  However,  a 
polarization  at  uj  —  5uj  also  results  due  to  the  beat  note  at  5uj  which  modulates  the  intensity: 
the  resulting  modulation  of  the  population  inversion  with  time  (i.e.,  a  population  pulsation) 
leads  to  nonlinear  frequency  mixing.  We  express  the  full  polarization  as 


p(t)  =  + c-c- 


(26) 

m=— ,0,+ 

where  c =  uj  +  Suj  and  cj_  =  uj  —  Suj.  We  can  solve  for  the  polarization  as  done  in  [3], 
keeping  only  terms  to  first  order  in  the  weak  field  £+,  which  gives 

£o 


_  ie 

^0  =  - OLWq- 

UJba  1 


^4 


^  16 
V+  =  — aw0 

^ba 

if  -  ~ 

V-  =  — aw0AtS0S0S: 

Uba 


[1  -i(A-Suj)T2 


A^o  S*S, 


(27) 

(28) 
(29) 


where 


At  = 


—  [1  +  i(A  +  fiw)T2](l  +  zAco’2'2  /  2) 


[1  -  i( A  -  Su)T2](1  +  iAT2)  [(1  +  iSuTJl  1  +  i( A  +  Su)T2][  1  -  i( A  -  Scu)T2]  +  (1  +  iSuT2) \S0\2 

(30) 


A+  = 


-(1  -  iSuT2/2) 


(1  -  iAT2)  (1  -  i8uTx)[  1  +  i{ A  -  5oj)T2}[  1  -  i( A  +  Su)T2]  +  (1  -  i5cuT2)\S0\2 


(31) 


are  the  self-mixing  and  cross-mixing  coupling  coefficients,  respectively  We  consider  the 
dipoles  to  be  embedded  in  a  host  medium  of  permittivity  e  and  permeability  fi.  (We  adopt 
the  convention  of  [4]:  e,  //  and  the  speed  of  light  c  =  l/  W/l  always  take  their  values  in 
the  background  host  medium.)  Many  of  the  material  properties  of  the  two- level  system  are 
lumped  into  the  “Beer  loss  rate” 

N  d2T2ojbac^/jr/e 


a  = 


h 


(32) 


which  is  related  to  the  more  familiar  Beer  absorption  coefficient  a  (with  units  of  inverse 
length)  that  appears  in  Beers  law  of  absorption  by  a  =  ac.  (Note,  however,  that  in  our 
expressions  for  the  polarization  due  to  the  two-level  system,  all  factors  of  e  and  /i  drop 
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out;  that  is,  these  expressions  do  not  contain  the  polarization  contributions  due  to  the 
background  medium.)  The  central  mode  amplitude  £0  has  been  normalized  such  that  £0  = 
/vv/TiT2£0.  Note  that  Vo  is  unaffected  to  first  order  in  £+.  The  polarization  V+  comes  from 
two  contributions.  First,  there  is  the  linear  contribution  from  the  Lorentz  oscillator  which 
£+  would  induce  even  in  the  absence  of  the  strong  field  £o .  Second,  there  is  a  contribution 
due  to  the  PP  which  is  described  by  the  term  A+.  The  term  V-  is  due  solely  to  the  PP  and 
is  governed  by  At.  Note  that  the  full  polarization  is  directly  proportional  to  the  steady- 
state  population  inversion  w0]  this  will  be  important  when  we  generalize  our  results  to 
standing-wave  cavities,  where  w0  varies  with  position. 

Now  that  we  have  the  polarization,  we  can  calculate  the  gain  seen  by  the  sideband  field. 
We  define  the  gain  g  (with  dimension  of  frequency)  of  the  sideband  as  the  power  density 
generated  at  u  +  6u>  by  the  interaction  of  the  field  with  the  dipoles-considering  only  field 
and  polarization  terms  oscillating  at  uj  +  ^-divided  by  the  energy  density  of  the  exciting 
sideband  field,  or 


_(EP)± 

2£|£+|2 

iu+(£+Vt 


£*+V+) 


2e|T+|2 


(33) 

(34) 


1.  A  =  0 


Here  we  consider  the  case  of  zero  detuning,  A  =  0,  which  simplifies  the  mathematical 
expressions  considerably  and  is  a  prerequisite  to  understanding  the  case  of  non-zero  detuning. 
Under  this  simplified  scenario,  we  denote  the  self-mixing  coefficient  by  A,  where 
^  _  —  (1  +  iSu>T2/2) 

(1  T  i8ujT\){\  T  idujTt)^  T  (1  +  i&nT^ITol2 

and  it  is  simple  to  show  that  the  cross-coupling  coefficient  At  is  simply  A*.  The  gain  of  the 
sideband  field  is  found  to  be 


g+  =  aw0 


irpf+w<A|141 


(36) 


(We  have  used  (cu+5cu)/u)ba  «1.)  Thus,  the  gain  can  be  nicely  divided  up  into  a  contribution 
from  the  Lorentz  oscillator  and  a  contribution  from  the  PP.  All  of  this  is  proportional  to 
aw 0:  a  gives  you  the  gain  of  a  weak  field  tuned  to  line-center  in  a  perfectly  inverted  medium 
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(or  alternatively,  the  loss  seen  by  a  weak  field  tuned  to  line-center  in  a  material  in  its  ground 
state),  and  w0  gives  you  the  expectation  value  of  finding  an  electron  in  the  excited  state 
(equal  to  1  when  excited,  -1  when  in  the  ground  state,  and  0  at  transparency).  Note  that 
Real  (A)  can  be  positive  or  negative,  which  we  will  discuss  shortly. 


B.  Three- frequency  operation 


Of  course,  the  polarization  created  at  u>  —  5cu  will  create  a  field  at  that  frequency,  which 
is  precisely  why  in  the  experiments  we  always  observe  the  two  sidebands  appearing  simul¬ 
taneously.  One  sideband  cannot  exist  in  isolation  when  the  mixing  terms  naturally  couple 
them  together.  Therefore,  we  need  to  consider  the  field 

£  =  £ o  +  £+ei5ut  +  £_e~iSi0t.  (37) 


The  polarization  at  each  sideband  frequency  now  contains  a  Lorentzian  term,  a  self-mixing 
term,  and  a  cross-mixing  term: 


_  xe 

V+  =  — aw0 

M ba 

—  ie  _ 

V _  =  — awo 

^ ba 


£, 


[1  -  i(  A  -  Slo)T2 ] 

£_ 

[1  -  i( A  +  Slj) T2] 


KX  £)£()£+  +  A+T0^o^- 

\-j0£*0£-  +  aiw; 


(38) 

(39) 


where  A_  and  A+  are  obtained  by  making  the  substitution  Su  —Soj  in  the  expressions  for 
A+  and  Al,  respectively,  given  in  Eqs.  30-31  . 


1.  A  =  0 


Let  us  again  focus  on  the  case  A  =  0,  for  which  the  polarization  at  each  sideband 
simplifies  to 


V+  =  — aw0 

^ ba 


V-  =  — awo 

Uba 


|t  +  V £»£'„£-  +  £1 


(40) 

(41) 


where  A  is  simply  Aljl  evaluated  for  A  =  0.  We  see  the  nice  property  that  when  A  =  0, 
Aljl  =  A+  (=  A),  and  A I  =  At  (=  A*);  in  other  words,  the  self-  and  cross- mixing  coupling 
coefficients  are  equal. 
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The  gain  g+  of  the  positive  sideband  is 


g+  =  aw0 


1 

i  +  {SuT2y 


+  Real 


(42) 


and  a  similar  expression  holds  for  the  minus  sideband.  This  equation  tells  us  that  the  PP 
contribution  to  the  gain  depends  on  the  phase  and  amplitude  relationships  of  So,  S-,  and 
S. |_,  which  is  not  too  surprising  because  the  amplitude  of  the  PP  itself  is  sensitive  to  these 
parameters.  Without  loss  of  generality,  we  can  take  So  to  be  real.  If  S+  =  £*,  then  the 
two  sidebands’  contributions  to  the  beat  note  at  5u>  add  constructively,  resulting  in  a  field 
whose  amplitude  modulation  (AM)  is  twice  the  strength  of  a  field  with  only  one  sideband. 
If  S. |_  =  —  S*_,  then  the  two  sidebands’  contributions  to  the  beat  note  at  Su  destructively 
cancel  and  there  is  no  longer  any  amplitude  modulation  at  frequency  Suj .  We  refer  to  such 
a  field  as  frequency- modulated  (FM).  We  see  from  Eq.  42  that  the  AM  sidebands  therefore 
experience  a  PP  contribution  to  the  gain  that  is  twice  as  large  as  the  single  sideband  case, 
while  the  FM  sidebands  experiences  only  the  background  Lorentzian  gain,  consistent  with 
the  fact  that  there  is  no  PP  in  this  case.  We  summarize  this  with  the  formula  for  the  gain 
g  of  each  sideband  for  the  case  of  equal- amplitude  sidebands  (|£+|  =  |£_|), 


g  =  aw0 


1  +  ((Lt2)2  +  Real(A)|£’« 


2  ;  AM 
0  ;  FM 


(43) 


Note  that  for  a  superposition  of  AM  and  FM,  the  gain  due  to  the  PP  will  fall  between  0 
and  2  times  the  factor  Real(A)|£0|2- 


VI.  THEORY:  INSTABILITY  THRESHOLD 

This  section  gives  a  more  detailed  derivation  of  the  instability  threshold,  combining  the 
populating  grating  and  the  population  pulsations  into. 

When  a  continuous- wave  (cw)  laser  is  pumped  at  its  lasing  threshold,  only  a  single 
frequency  of  light-the  one  nearest  the  gain  peak  that  also  satisfies  the  roundtrip  phase 
condition-has  sufficient  gain  to  overcome  the  roundtrip  loss  and  begins  to  lase.  As  the 
pumping  is  increased,  the  single-mode  solution  yields  to  multimode  operation;  this  is  known 
as  the  single- mode  instability.  Our  goal  is  to  determine  1)  how  hard  to  pump  the  laser  to 
reach  the  single-mode  instability  and  2)  which  new  frequencies  start  lasing. 
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Consider  a  laser  pumped  above  threshold  that  is  lasing  on  a  single-mode,  which  we 
refer  to  as  the  primary  or  central  mode.  If  another  mode  is  to  lase,  it  must  be  seeded  by 
a  spontaneously  generated  photon  at  a  different  frequency.  This  photon  will  necessarily 
create  a  beat  note  through  its  coexistence  with  the  primary  mode,  resulting  in  a  population 
pulsation.  The  gain  seen  by  the  new  frequency  must  therefore  account  for  this  parametric 
gain  in  addition  to  the  background  Lorentzian  gain.  Furthermore,  the  PP  couples  the 
sideband  to  the  symmetrically  detuned  sideband  frequency  on  the  other  side  of  the  primary 
mode,  so  we  should  in  general  assume  the  presence  of  both  sidebands.  Because  the  instability 
threshold  depends  on  the  cavity  geometry,  we  will  consider  a  traveling-wave  laser  as  well 
as  a  standing-wave  laser.  In  both  cases,  the  strategy  is  the  same.  First,  we  solve  for  the 
single-mode  intensity  £q  and  the  population  inversion  w$ (z)  as  a  function  of  the  pumping, 
entirely  neglecting  the  sidebands.  Knowing  this,  we  can  then  calculate  the  sideband  gain  in 
the  presence  of  the  primary  mode. 

We  start  with  the  wave  equation 

d2E  1  d2E  _  d2P 

Following  the  approach  used  to  calculate  the  optical  parametric  oscillation  threshold  in 
optically  pumped  microresonators  [5],  we  expand  the  field  in  terms  of  the  cold  cavity  modes, 

E{z,t)  =  ^2  £m(t)T  m{z)elUmt  +  C.C.  (45) 

m=— ,0,+ 

The  spatial  modes  obey  the  normalization  condition 

J  dz  \Tm(z)\2  =  1.  (46) 

When  group  velocity  dispersion  (GVD)  is  non-zero,  the  two  modes  cu+  and  c 0-  will  not  be 
equidistant  from  ay.  We  have  also  assumed  that  the  spatial  and  temporal  dependence  of  the 
modes  can  be  separated.  This  is  a  good  approximation  in  the  case  of  a  laser,  because  we 
know  the  intracavity  field  will  be  sharply  resonant  at  the  modes.  The  spatial  variation  of  the 
polarization  can  be  described  by  making  the  substitution  £m  — >  £mTm(z)  and  w0  —>  w0(z) 
into  the  polarization  Eqs.  38-39,  which  results  in  the  polarization 

p(z,t)=  J2  'PmMe^  +  c.c.  (47) 

m=— ,0,+ 
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where 


-^—aw0(z) 

^ ba 

£+T  +(J) 

[1  -  i{A  -  5cu)T2] 

—  mi)nM 

T_T  _{z) 

Uc  (JJ  o  l  sO  ) 

Uba 

.[!-*(  A  +  Soo)T2] 

+  Ai(3)|T„(2)|2T+(2)|£„|2£+  +  A;(3)T(l(3)2T:(3)e”*42£: 

(48) 

+  A:(z)|T„(z)|2T_(z)|£0|2£_  +  Ai(s)T„(z)2T;(z)e“£2£; 


(49) 


We  have  introduced  u  =  2cj0  —  u>+  —  u>-,  the  deviation  of  the  cold  cavity  modes  from  equal 
spacing.  Note  that  the  As  now  depend  on  z  due  to  the  term  in  their  denominators  dependent 
on  the  primary  mode  amplitude.  Because  we  no  longer  demand  that  the  two  sidebands  have 
the  same  detuning  Suj ,  and  At  should,  strictly  speaking,  be  calculated  using  the  detuning 
Suj. |_  =  uj+  —  wo,  while  A I  and  A+  should  depend  on  Scu-  =  ujo  —  tu_.  In  practice,  we  can 
ignore  this  difference  in  the  As;  the  term  elujt  captures  the  most  important  effect  of  GVD. 

Plugging  everything  into  the  wave  equation  gives 


m 


A  ^2  -  /Pm,]OBB)etUmt 

m 

(50) 


where  the  slowly-varying-envelope  approximation  allowed  us  to  ignore  second  time  deriva¬ 
tives  of  Sm  on  the  left-hand  side,  and  first  and  second  derivatives  of  £m  on  the  right-hand 
side.  The  spatial  modes  T m(z)  are  chosen  so  that  the  first  term  on  the  LHS  equals  zero.  The 
loss  of  each  mode  has  been  added  to  the  equation  in  the  form  of  a  polarization  contribution; 
we  assume  each  mode  has  the  same  linear  loss,  which  can  be  expressed 


Z€  — 

Vm,\oss{z,t)  =  - £T  m(z)£m(t). 

^ba 


(51) 


Equation  50  couples  all  of  the  modes  £m.  We  can  project  this  equation  onto  each  mode 
by  multiplying  by  T n(z)  and  integrating  over  the  length  of  the  laser  cavity,  thus  taking 
advantage  of  the  orthonormality  of  the  spatial  modes  T m(z),  and  then  equating  terms  which 
oscillate  at  the  same  frequency  (since  terms  with  different  frequencies  will  not  affect  the 
time-averaged  gain  seen  by  a  mode).  The  result  is  one  equation  for  the  central  mode 


a 

2(1  -  iATt) 


h  «,„(*)|t„mi2 


(52) 
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one  for  the  positive  sideband 


•  £  a 

£+  =  ~2€++2 


£+ 


l  -  i(A  -  5lu)T2  J  L 


H  7 

-  Wq(z)\T  +(z)\2 


+  \£o\2£^ 


dz 

T 


w0(z)AX(z)\r0(z)\2\r+(z)f 


dz 


+£2£le^  /  —  wo(z)A-(z)To(zyT*_(z)T*+(z) 


and  one  for  the  negative  sideband 


•  ^  „  a 

£-  =  -2£-+2 


£4 


T-  wo(v>|T_(V)|2 


1  —  i(A  +  5uj)T2  J  L 
r  dz 

7\£o?£-  -  w0(z)AZ(z)\T0(z)\2\T.(z)\^ 


L 

dz 


+£l£*+eiM  /  —  w0(^)A1(^)T0(^)2T;(^)T1(^) 


(53) 


(54) 


These  three  equations  will  be  used  to  understand  the  instability  threshold.  In  general,  one 
must  first  apply  the  steady-state  condition  £0  =  0  to  Eq.  52  which,  together  with  the  Bloch 
equation  relating  the  field  to  the  inversion,  will  yield  the  amplitude  of  the  primary  mode  £0 
along  with  the  resulting  population  inversion  w0(z),  both  as  a  function  of  the  pumping  weq. 
This  information  is  then  used  in  Eqs.  53-54  to  determine  the  minimum  level  of  pumping  weq 
at  which  a  pair  of  sidebands  with  detuning  Su  experiences  more  gain  than  loss.  This  is  the 
instability  threshold. 

So  far,  we  have  kept  Eqs.  52-54  as  general  as  possible  to  account  for  arbitrary  spatial 
profiles,  GVD,  and  detuning  A  between  the  lasing  mode  and  the  peak  of  the  gain  spectrum. 
From  here  on  we  will  simplify  the  problem  by  taking  Q  =  0  (zero  GVD)  and  A  =  0,  and 
apply  these  conditions  to  the  simplest  possible  traveling-wave  and  standing-wave  cavities. 


A.  Traveling- wave  cavity 

For  the  traveling- wave  laser,  the  spatial  modes  are 

T  m(z)  =  e~ikmZ  (55) 

so  every  point  in  the  cavity  sees  the  same  intensity.  At  and  above  threshold,  the  population 
inversion  is  everywhere  saturated  to  the  threshold  inversion,  so  w0  is  independent  of  z.  For 
A  =  0,  the  inversion  is 

=  wth  =  rr  (56) 

a 
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and  the  intensity  of  the  primary  mode  is  given  by 


\£0\2=P-i 


(57) 


where  we  have  made  use  of  the  normalized  primary  mode  amplitude  So  =  h;\/T\T2So.  and  p 
is  the  pumping  parameter  defined  as  p  =  weq/wth.  Because  So \ 2  is  independent  of  z,  all  of 
the  As  are  independent  of  z.  Furthermore,  since  both  w0  and  the  As  are  independent  of  z, 
they  can  be  pulled  out  of  the  spatial  integrals  in  Eqs.  53-54.  These  integrals  are  then  equal 
to  one,  where  we  have  used  the  zero  GVD  condition  Co  =  0  in  order  for  the  cross-overlap 
integral  (the  last  integral  in  each  equation)  to  equal  one.  The  sideband  equations  become 


£+  = 


S  —  --S 


awth 

2 

awth 


S+ 


1  T  iSujT2 

S- 

1  —  %5uoTo 


+  k\So\2S+  +  A  SlS*_ 

+  a*|t0|2t_  +  a*42t; 


which  can  be  written  in  matrix  form 


M+  R+\  S 4 
R*_  M*_  /  \  S*_ 


where 


M+  =  M*_ 

R+  =  R*_ 


~2  + 


awth 


ih 


A|£0| 


( — - — 
y  1  H-  i5cuT2 


(58) 

(59) 

(60) 

(61) 

(62) 


(In  the  last  step,  we  have  finally  taken  the  freedom  to  choose  So  to  be  real,  which  we  can  do 
at  this  point  without  loss  of  generality.) 

Now,  if  we  assume  a  solution  of  the  form  S±  ~  eAt,  we  find  the  two  solutions  for  A 


A  =  i [M+  +  Ml  ±  yj(M+  -  Mlf  *  4B+R1], 


(63) 


The  net  gain  seen  by  each  sideband  is  given  by  Real(2A)  (the  factor  of  two  is  for  intensity 
gain  rather  than  amplitude  gain),  which  includes  the  gain  minus  the  loss.  Subtracting  off 
the  loss,  the  gain  g  seen  by  each  sideband  is 


g  =  awth 


1 

1  +  (5o;r2)2  + 


2  Real(A)|T0|2  ;  AM 
0  ;  FM 


(64) 
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FIG.  S5.  The  sideband  gain  g  am / £  of  a  traveling-wave  laser,  given  in  Eq.  66,  is  plotted  at  various 
pump  strengths,  for  three  different  values  of  Z:  1,  10,  and  100.  The  largest  value  of  p  in  each  plot 
is  equal  to  the  instability  threshold  given  in  Eq.  69. 


where  the  two  solutions  correspond  to  AM  and  FM  sideband  configurations.  Finally,  we 
recognize  that  the  gain  is  pinned  at  threshold,  so  awth  =  I,  and  we  write  down  the  sideband 
gain  normalized  to  the  loss 


£ 

£ 


1  +  (<Lt2)2  +Real(A)|f»|2 


2  ;  AM 
0  ;  FM 


(65) 


When  the  gain  g  exceeds  the  loss  £,  the  weak  sideband  amplitudes  experience  exponential 
growth,  therefore  the  single-mode  solution  becomes  unstable.  Note  that  the  Lorentzian  term 
is  always  less  than  1.  This  is  a  direct  result  of  uniform  gain  clamping  in  the  traveling- 
wave  laser,  which  clamps  the  net  gain  of  the  mode  at  the  peak  of  the  Lorentzian  to  zero, 
and  therefore  any  mode  detuned  from  the  peak  will  see  slightly  more  loss  than  gain.  FM 
sidebands  therefore  never  become  unstable  because  they  only  see  the  Lorentzian  gain.  On 
the  other  hand,  AM  sidebands  induce  a  PP  and  with  it  a  coherent  gain  term,  which  can 
provide  enough  extra  gain  on  top  of  the  Lorentzian  background  to  allow  the  sidebands  to 
lase, 


qv  =  T+jtjMf  +  2  Real(A)|t>|2-  (“) 

To  get  a  feel  for  the  sideband  gain,  we  have  plotted  Qam/£  in  Fig.  S5  at  various  pump 
strengths  p  for  Z  =  1,  10,  and  100,  where  Z  =  Ti/T2.  Graphically,  we  see  that  at  large 
enough  p  sidebands  will  become  unstable.  Analytically,  it  is  a  simple  matter  to  calculate 
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how  hard  to  pump  the  laser  p  before  the  sidebands  appear,  starting  from  Eq.  66.  We  start 
by  replacing  |£0|2  with  p  —  1,  and  note  that  this  substitution  must  also  be  made  in  A, 
which  implicitly  varies  with  | £0 1 2 .  Then,  setting  Qa m / 7  equal  to  one,  we  can  solve  a  simple 
quadratic  formula  for  Slu2, 


(5luT2)2 


-1  +  3 Zip  -  1)  ±  y/[l  -  3 Zip  -  l)]2  -  8 ZMp  -  1) 

2Z2 


(67) 


Finally,  we  must  apply  some  physical  reasoning:  as  p  is  increased  past  1,  the  sideband  gain 
increases.  Right  at  the  moment  when  the  instability  threshold  is  reached,  Slu2  must  take 
on  a  single  value.  Thus,  we  set  the  radical  in  Eq.  67  to  zero  and  solve  for  p.  After  solving 
another  simple  quadratic  equation,  we  find  that 

P  =  5+f±q/l  +  Jt  +  ^.  (68) 

How  do  we  choose  between  the  plus  and  minus  sign?  By  plugging  this  expression  for  p  back 
into  Eq.  67,  it  is  simple  to  check  that  only  the  plus  sign  yields  real- valued  solutions  for  Slu. 
Thus,  we  have  found  the  instability  threshold,  which  we  denote  Prngh , 

Prngh  =  5  +  |  +  4\A+^  +  ^5  (69) 

because  it  is  the  well-known  instability  threshold  found  by  Risken  and  Nummedal  (see  Eq. 
3.10  in  [6])  and  Graham  and  Haken  (see  Eq.  7.35  in  [7]).  Plugging  this  value  of  p  into  Eq. 
67  yields  the  value  of  Sou  of  the  sidebands  when  the  instability  sets  in 

(SurnghT2)2  =  ±  +  |  ^1  +  ^1  +  J  •  (70) 

One  thing  to  notice  is  that  in  the  limit  Z  1  (transverse  relaxation  must  faster  than 
longitudinal  relaxation),  the  instability  threshold  Prngh  —>  9  from  above  and  SujrnghT2  — > 
^12/Z. 


B.  Standing- Wave  Cavity 

As  before,  we  restrict  ourselves  to  the  case  A  =  0  and  lu  =  0.  We  will  see  that  calculations 
for  the  standing-wave  cavity  are  significantly  more  complicated  than  for  the  traveling-wave 
cavity.  The  spatial  variation  of  the  primary  mode  causes  the  inversion  w0  and  the  coupling 
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A  to  both  depend  on  z,  which  makes  the  integrals  more  difficult  to  compute.  For  this  reason, 
we  treat  the  problem  to  first  order  in  the  primary  mode  intensity  |£0|2,  which  allows  us  to 
compute  the  integrals  analytically.  However,  the  theory  can  be  extended  to  higher  order  at 
will,  or  the  integrals  can  always  be  computed  numerically. 

For  the  standing-wave  laser  with  perfectly  reflecting  end  mirrors,  the  spatial  profile  of 
each  mode  is  given  by 

T  m(z)  =  V2COS  (kmz).  (71) 


The  spatial  modulation  of  the  intensity  is  responsible  for  the  spatial  modulation  of  the 
population  inversion  wo(z),  though  mitigated  somewhat  by  carrier  diffusion.  We  calculated 
wq(z)  in  Sec.  IV  of  the  supplementary  information.  The  result  is 


w0(z)  =  Wth 


1  + 


7  D  P~  1 
2  1  +  7d/2 


p  —  1 

~7dT+W2 


cos(2£r0£) 


(72) 


to  first  order  in  |£0|2,  where  wth  =  f/7,  and  7 d  =  (1  +  4k2DTi)~1  is  the  diffusion  parameter. 
The  spatial  variation  of  the  inversion  has  important  consequences.  For  one,  it  reduces  the 
power  of  the  laser,  which  is  given  by 


p  —  1 

1  +  7i>/2 

Secondly,  the  gain  is  no  longer  uniformly  clamped  by  the  primary  lasing  mode,  which  will 
allow  new  modes  to  lase  even  in  the  absence  of  PPs. 

The  spatial  variation  of  the  primary  lasing  mode  also  causes  A  to  vary  with  position.  In 
keeping  with  our  approximations,  we  can  expand  A  to  zeroth  order  in  |£0|2  because  in  our 
equations  A  always  multiplies  |£0|2,  so  the  final  result  is  first  order  in  \£q\2 .  We  define  the 
zeroth  order  expansion  of  A  to  be 


(3)  — (1  +  hW'77/2) 

X  =  (l  +  hkjTi)(l  +  hkjT2)2’ 

where  the  symbol  y(:b  was  chosen  to  emphasize  that  this  term  now  plays  the  role  of  a 
third-order  nonlinear  coefficient. 

We  start  with  the  sideband  Eqs.  53-54,  replace  wo(z)  with  Eq.  18,  A (z)  with  and 
keep  only  terms  to  first  order  in  |£0|2.  The  resulting  equation  for  the  growth  of  the  positive 
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sideband  is 


f+  =  —  A4 


i  +  yifoi2 

1  -|-  iSujT2 


+  x{3)\So\2£+J  ^  |ToW|2|th 


+x(3)42£- 


and  a  similar  equation  can  be  written  down  for  £_.  We  define  the  longitudinal  overlap 
integrals 


rsfiif  — 


\Tn(z)\2\Tjz)\2  =  1 


r 

-1-  cross 


rjy 

l  TM^mz)T*+(z)  =  1/2. 


The  implication  is  that  the  self-mixing  interaction  of  a  sideband  with  itself,  mediated  by 
the  primary  mode  intensity,  is  twice  as  large  as  the  cross-mixing  interaction  of  one  sideband 
generating  gain  for  the  other  sideband,  again  mediated  by  the  primary  mode  intensity.  This 
is  true  only  for  the  cosine-shaped  modes  that  we  have  assumed,  and  the  overlap  integrals 
will  change  when  the  longitudinal  spatial  profile  changes,  as  when  the  non- unity  reflectivity 
of  the  facets  is  taken  into  account.  The  sideband  Eqs.  53-54  become 


OLWth  1  +  ^|£o|2 

2  IT  iSuT2 


£+  +  rselfX(3)|£0|2£+  +  rcrossX(3)42r 


awth  h  +  ^i^oi2 


1  —  iScuT \ 


+  rsel{X(3r\£o\2£-  +  rcrossx(3)*42^ 


which  we  express  as 


where 


M+  =  M*  =  — 
+  “  2 


M+  R+ 
R*  M* 


£  OWth  (  1  T  '^2'  | Sq  | 2 
2  +  2  l  1  T  iSojT2 


rseifX(3¥o|2 


R+  =  R*_  =  ^(TCIOSsX{3)\£o\2)- 


As  we  did  for  the  traveling-wave  laser,  the  sideband  gain  is  easily  calculated  from  these  two 
coupled  first-order  differential  equations.  Normalizing  the  gain  to  the  total  loss,  we  find 


g  IT 


£  1  T  (5uT2)* 


9  +  Real[x(3)]|<£o 


Ueif  +  rcross  —  |  ;  AM 

r seif  —  rcroSs  =  \  ;  FM 
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There  are  two  things  to  notice  here.  As  the  laser  pumping  is  increased,  the  term  7d|£o|2/2 
grows,  and  consequently  the  gain  is  not  clamped  at  the  threshold  value.  This  is  due  to  spatial 
hole  burning,  or  more  precisely,  the  imperfect  overlap  of  the  standing-wave  modes  together 
with  a  finite  amount  of  carrier  diffusion.  We  view  this  background  gain  as  a  Lorentzian- 
shape  whose  amplitude  increases  with  the  pumping,  and  is  therefore  fully  capable  of  pulling 
the  sidebands  above  threshold,  without  any  additional  PP  contribution  to  the  gain. 

Secondly,  the  PP  contribution  to  the  gain  never  vanishes.  Even  when  the  sidebands  are 
phased  such  that  an  FM  waveform  is  emitted  from  the  laser,  there  is  still  a  PP  within  the 
laser  cavity.  The  reason  for  this  is  the  imperfect  overlap  of  the  two  sidebands’  spatial  modes, 
which  means  that  at  any  given  position  within  the  cavity,  the  plus  and  minus  sideband  are 
likely  to  have  different  amplitudes.  Therefore,  even  if  the  two  sidebands  are  phased  such 
that  their  contributions  to  the  beat  note  at  8u  destructively  interfere  with  each  other,  the 
destruction  is  not  perfect.  The  amplitude  of  the  PP  varies  with  position  in  the  cavity,  and 
in  locations  where  the  two  sideband  amplitudes  are  equal  the  PP  will  not  exist,  but  the 
spatially  averaged  effect  of  the  FM  PP  yields  the  factor  of  1/2  in  Eq.  83.  By  the  same 
token,  sidebands  phased  for  AM  will  not  fully  constructively  interfere,  yielding  a  factor  of 
3/2  for  the  PP  contribution  to  the  gain  rather  than  the  factor  2,  as  it  would  be  for  the 
traveling- wave  laser. 


VII.  NUMERICALLY  CALCULATING  THE  INSTABILITY  THRESHOLD 

In  this  section,  we  demonstrate  the  predictions  of  the  theory  for  the  three  uncoated  lasers, 
and  compare  the  results  with  the  measurements. 

Because  the  theory  assumes  end  mirrors  with  unity  reflectivity,  we  can  only  expect  Eq. 
83  to  apply  reasonably  well  to  the  uncoated  QCLs.  For  each  device,  ~/d  is  calculated  using 
the  theoretical  value  of  Tup  (calculated  from  the  bandstructure)  and  the  diffusion  constant 
D  =  77  cm2/s  [8],  giving  7d  =  0.4  (DS-3.8),  0.49  (TL-4.6),  and  0.93  (LL-9.8).  For  these 
large  values  of  7d,  the  incoherent  gain  increases  rapidly  with  the  pumping,  and  we  find  from 
Eq.  83  that  the  FM  instability  will  have  a  lower  threshold  than  the  AM  instability,  regardless 
of  the  value  of  Tj.  The  gain  recovery  time  Tj  of  each  QCL  is  not  as  easily  calculable  as  Tup 
because  it  depends  on  a  few  other  time  constants  of  the  active  region,  such  as  the  escape  time 
of  the  electron  from  one  injector  region  to  the  next  active  region.  Therefore,  we  treat  Tj  as  a 
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Ti/T2 


FIG.  S6.  Numerical  solutions  of  the  instability  threshold  obtained  by  setting  the  gain  g  in  Eq.  83 
equal  to  the  loss  i,  yielding  both  (a)  the  sideband  separation  and  (b)  the  pumping  psi:).  The 

experimentally  measured  values  of  are  compared  to  the  theory  to  infer  T1/T2,  which  also 

gives  the  theoretical  prediction  for  the  instability  threshold  ps\,. 


variable  and  calculate  the  instability  threshold  ps b  and  sideband  spacing  5cysb  as  a  function  of 
T\.  The  resulting  curves  are  shown  in  Fig.  S6.  By  comparing  the  curves  with  the  measured 
values  of  we  can  deduce  the  values  Ti  =  1.83  ps  (DS-3.8),  1.15  ps  (TL-4.6),  and  0.91 
ps  (LL-9.8).  For  these  values  of  Ti,  the  theory  predicts  an  instability  threshold  of  ps b  =  1.02 
(DS-3.8),  1.09  (TL-4.6),  and  1.04  (LL-9.8).  It  is  encouraging  that  these  fitted  values  of  Tf 
are  close  to  the  accepted  value  of  the  QCL  gain  recovery  time,  which  has  been  shown  by 
pump- probe  experiments  [9,  10]  and  theory  [11]  to  be  around  2  ps.  However,  the  predicted 
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psb  is  significantly  lower  than  the  measured  values  Jsb/-Jth  =  1.12  (DS-3.8),  1.17  (TL-4.6), 
and  1.14  (LL-9.8),  and  the  discrepancy  is  made  worse  by  the  fact  that  J/Jth  is  likely  an 
underestimate  of  p  (see  the  discussion  in  Sec.  I  of  the  supplement).  The  fact  that  the  theory 
underestimates  the  instability  threshold  is  perhaps  not  surprising,  as  we  have  only  made 
sure  that  one  of  the  two  necessary  conditions  for  sideband  oscillation  is  satisfied  (gain,  not 
phase).  We  hope  that  future  work  which  accounts  for  the  detuning  A,  the  detuning  between 
the  lasing  mode  and  the  cold  cavity  mode  it  occupies,  and  GVD  can  accurately  predict  the 
instability  threshold,  which  would  be  a  milestone  in  the  understanding  of  lasers,  and  also 
yield  a  novel  laser  characterization  method  of  lifetimes  and  diffusion  rates  by  comparing 
measured  values  of  ps b  and  5us b  to  an  established  theory. 
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